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ABSTRACT 

Dust and associated depletion of heavy elements from the gas phase can modify the 
thermal properties of H II regions from the dust-free case, with significant consequences 
for the emergent optical spectrum. We present the results of theoretical calculations 
illustrating the effects of grains on the spectra of giant, extragalactic HII regions, with 
emphasis on high metallicity systems (Z £ Z ). Dust provides a simple explanation for 
the observational absence of pure Balmer-line spectra that are expected on theoretical 
grounds for dust-free, chemically enriched nebulae. Grains may also play a role in 
enhancements of forbidden- line emission observed in HII regions in the enriched nuclei 
of normal galaxies. In most cases, depletion introduces the strongest perturbations 
to the optical spectrum. Selective absorption of the ionizing continuum as well as 
heating by grain photoelectrons are important in some instances, however, and grain 
heating can be particularly important for enhancing emission in high- ionization lines. 
Allowing for depletion, the presence of dust is unlikely to introduce large errors in 
global metallicity indicators, although uncertainties in depletion factors coupled with 
the sensitivity of infrared cooling to electron density will make accurate calibrations 
difficult at high Z. The present calculations establish further that previous relative 
abundance analyses that fail to take into account dust effects in a self-consistent way 
(grain heating as well as depletion) may overestimate temperature gradients in high-Z 
nebulae, resulting in errors in relative abundances for different elements. 



Subject headings: ISM: H II Regions - ISM: Dust - ISM: Abundances 



1. Introduction 



Giant H II regions have a fundamental role as probes of elemental abundances in extragalactic 
environments. Existing spectroscopic studies of extragalactic HII regions tend to emphasize 
systems with subsolar metallicities, due in part to inherent interest in chemically unevolved 
systems, but also because of selection effects, since characteristic luminosities and emission-line 



bubble Fellow 
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equivalent widths for HII regions tend to decrease with increasing source metallicity (e.g., 



Kennicutt 1988). Data for nebulae in relatively high- metallicity environments, such as the inner 
regions of giant galaxies, are becoming increasingly available, however, and additional theoretical 
work is necessary if we are to make use of this information for quantitative abundance analyses. 

The presence of abundant atomic coolants in enriched HII regions will typically lead to a 
low equilibrium electron temperature, with a large fraction of the radiated energy emitted by 
the heavy elements through infrared ground-state fine structure transitions. The strengths of the 
comparatively energetic optical lines are very sensitive to the electron temperature, and hence 
to nebular parameters that affect temperature. At abundances above a few times solar, the 
fine-structure cooling is efficient enough that optical forbidden lines from simple, low-density H II 
regions are expected to be unobservable (e.g., Rubirj 1985). To date, however, no HII region 



emitting a pure Balmer-line spectrum has been discovered, even in surveys of Ha-selected HII 
regions in environments for which metallicity Z ^> Z Q is expected ( Oey fc Kennicutt] 1993). The 



lack of such objects implies that physical parameters in real nebulae are modified so that the 
electron temperature is high enough to produce significant optical forbidden-line emission. 

The effects of interstellar dust could conceivably underlie at least some of the enhancement of 
optical emission in high-Z HII regions. Dust can act to modify the thermal structure of a nebula 
in several ways. The removal of atomic coolants from the gas phase via depletion will tend to 
increase the electron temperature, leading to measurable changes in the emitted spectrum (e.g., 
Henryj 1993). Grains may selectively absorb a fraction of the Lyman continuum photons powering 



the nebula, resulting in modification to the ionizing spectral energy distribution seen by the gas 
and to the resultant heating (e.g., Sarazin 1977; |Mathisj 1986; [Aannestad| 1989). Grains will also 



participate directly in the thermal balance by contributing photoelectrons that heat the plasma, 



and by cooling through radiation of energy transferred in captures of electrons (e.g., Oliveira & 



Maciel 1986; Baldwin et al. 1991). 



Previous studies emphasizing nebulae at Z £ Z Q have generally concluded that selective 
absorption of the ionizing continuum by dust has a relatively small effect on the emergent optical 
spectrum of diffuse HII regions (e.g., |Mathis| 1986). Given the greater temperature sensitivity of 
optical lines at higher Z, we might expect the influence of dust to be more important in metal-rich 
nebulae; this statement is especially true if the dust-to-gas ratio scales with Z. The observational 
consequences of grain heating and cooling within H II regions over the entire span of metallicity 
have also received only limited attention in previous work. 

Motivated by these considerations, we have undertaken photoionization calculations 
appropriate for giant extragalactic HII regions in order to understand possible consequences 
of dust in chemically enriched nebulae. We find that dust and associated depletion can indeed 
produce nonnegligible perturbations in the optical spectrum of metal-rich HII regions. The 
influence of dust in such objects provides a possible basis for explaining several unexpected 
empirical characteristics of HII regions in high-Z environments, including the lack of pure 
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Balmer-line nebulae. If grains survive in the inner parts of spherical nebulae, they may selectively 
enhance emission in high- ionization lines. Dust is unlikely to contribute large uncertainties in gas 
phase abundance estimates based on the oxygen emission lines, but ambiguities in depletion factors 
and density effects will contribute substantial uncertainty to "empirical" metallicity-excitation 
calibrations for Z ^ Z© (e.g., Pagel et~aL| 1979). Finally, we suggest possible revisions to 
theoretical relations between ion-weighted temperatures employed in relative abundance analyses, 
necessitated by a self-consistent treatment of the presence of dust in HII regions. 



2. Nebular Calculations with Dust 



Calculations appropriate for prediction of H II region optical spectra were generated with 



the photoionization code Cloudy, version 84.09 ( Ferland 1994). Parameters were chosen such 



that the results would illustrate the possible consequences of dust in giant, extragalactic HII 



regions (GEHRs) as a function of metallicity. The ionizing continuum was represented by Kurucz 
(1979) stellar atmospheres corresponding to T* of 38,000 K and 45,000 K, with total ionizing 
luminosities (photons s -1 ) of log Q(H) = 50.7 and 51.2, respectively. The cooler stellar continuum 
model is interpolated from theoretical calculations for T* = 35,000 K and 40,000 K, and both of 
the final continuum choices assume stellar log(g) = 4.5. Based on systematic trends in existing 
observations, the lower T* values are more applicable for nebulae with relatively high Z (e.g., 



McCall, Rybski, & Shields 1985; [Vflchez & Pagel| 1988). 



The nebular geometry was assumed to be spherical, extending from an inner radius of 0.01 
pc to the Stromgren radius R s defined as the point where the ionized hydrogen fraction fell below 
1%. A total hydrogen density of nn = 10 cm -3 was adopted for material described by volume 
filling factors of e = 0.0617 and 0.0347 for the 38,000 K and 45,000 K cases, respectively. If we 
define ionization parameter U in the customary way [U = Q (H)/(47r R 2 n e c), where c is the speed 
of light], these combinations of parameters yield log U = —3.0 at Z = Z in the dust-free case; 



estimates of U for GEHRs are comparable to or somewhat larger than this value (Shields 1990 
and references therein). U is a function of Z, due to the dependence of recombination rates and 
hence R s on electron temperature T e , as well as a function of dust content, since grains will absorb 
a portion of the ionizing radiation and also modify the T e structure of the nebula through heating 
and cooling processes. 



Total elemental abundances for the case of Z = Z & were taken from Grevesse fc Anders 
(1989). For other values of Z, abundances of elements heavier than helium were assumed to be 
in solar proportions, with the exception of nitrogen, which was assumed to have a secondary 
component such that its abundance scales as Z 2 . Empirical evidence for a secondary contribution 
to N in GEHRs is discussed, for example, in Vila-Costas fc Edmunds- (1993). Depletion onto grains 



was assumed to modify gas-phase abundances by factors taken from Gowie fc Songaila (1986, 
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Table 3) corresponding to an overall average for the interstellar medium (ISM). The resulting total 
and depleted abundances at Z = Z Q are listed in Table 1. 

The grains themselves were assumed to consist of a graphite-silicate mix with a size 



distribution as described by Mathis, Rumpl, fc Nordsieck (1977) for the diffuse ISM, with optical 
constants as given by Draine & Lee (1984) and extended within the extreme ultraviolet (EUV) 
bandpass by Martin &: Rouleau] (1989). The resulting dust opacity features a maximum at photon 
energies of ~ 17 eV and a decrease at higher energies, as required by the Kramers-Kronig relations 
for realistic grain materials (e.g., Bohren Sz Huffman 1983; Mathis- 1986). 



Scattering by grains is neglected in the present calculations. At ionizing energies, the 
scattering phase function for the grains is expected to be predominantly forward-throwing. At 
optical wavelengths, scattering in and out of the line of sight will approximately cancel in a 
spherical geometry. In the current implementation of Cloudy, reddening of the emergent spectrum 
by grain absorption is also neglected in the spherical case. Our principal interest in the present 
study is the consequences of dust for the structural and thermal characteristics of nebulae, rather 
than on emission-line transfer effects of the dust. The effects of internal scattering and absorption 
have been discussed previously by Mathisj (1970, 1983), |Petrosian &: Dana] (1980), and Aannestad 
(1989), among others. For the cases of interest here, the scattering and absorption optical depths 
within R s are relatively small (§ (4.3.|) ; in such cases, dereddening of the observed optical spectrum 
based on the Balmer line strengths is expected to reproduce the intrinsic line ratios with reasonable 
accuracy ( |Mathis| 1983). As a result, the predictions of line ratios presented here are suitable for 
direct comparison with dereddened ratios derived from observations. 

Thermal effects of dust are quantified by extending the calculation of nebular ionization and 
energy equilibrium to include grains. Additional details concerning the treatment of grain physics 
within Cloudy can be found in Baldwin et alj (1991, Appendix C) and Fcrland| (1994). The 
dust-to-gas ratio was assumed to scale in proportion with Z, with a normalization of 0.01 by mass 
at Z = Zq. 

Giant H II regions exhibit supersonic line widths that correlate in amplitude with emission-line 
luminosity ( Melnick] 1979). Line widths in excess of thermal values will modify the transfer of 
lines including infrared fine-structure transitions that would otherwise be optically thick in a 
quiescent nebula. The added velocity width facilitates the escape of these lines, and thus may 
decrease the equilibrium T e when the infrared transitions contribute a significant fraction of the 



total cooling (e.g., Simpson 1975). We consequently assumed the presence of a microturbulent 
velocity component bt ur b added in quadrature with the thermal b (= [2kT e /m] ' 5 ), with values of 
hurb = 17.4 km s -1 for the case of T* = 38, 000 K and log Q(R) = 50.7, and b turb = 23.0 km s" 1 for 
T* = 45,000 K and log Q(H) = 51.2. These choices of b tur b are derived from the empirical relation 
between recombination line luminosity and velocity width given by |Arsenault &TRxyy\ (1988). 
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3. General Results 



The predicted strengths of optical forbidden lines, normalized to hydrogen recombination 
features, are plotted as a function of Z in Figures 1 and 2 for T* = 38, 000 K and 45,000 K, 
respectively. The curves show the expected behavior with increasing Z, with forbidden-line 
emission initially rising in response to the growing number of emitters, reaching a maximum, 
and declining at high Z as the total cooling efficiency increases and T e diminishes. The ratio 
of [Oil] A3727/H/3 exhibits an additional reversal, with increasing values at very high Z. This 
behavior stems from photoionization of O into excited states of + that can subsequently decay 
radiatively. The resulting emission component consequently scales directly with the abundance 
of O and hence with Z. While such emission can be important in some contexts (e.g., novae; 
Ferland fc Truran] 1981), it is unlikely to be significant in HII regions, for which [OIl]/H/3 ^ 0.3 



observationally (e.g., Oey fc Kennicuti] 1993). 



The separate curves in a given plot in Figures 1 and 2 illustrate the change in emission 
line behavior as the effects of dust are added incrementally. The solid curve corresponds to 
dust-free, undepleted nebulosity, while the other curves show the results of introducing depletion, 
dust opacity, and grain heating + cooling. In almost all cases, the effects of dust lead to an 
enhancement of optical forbidden-line emission over the undepleted, dust-free case, due to an 
increase in equilibrium T e . 

For most of the lines, the largest shift in the plotted line ratio curves results when depletion 
is introduced. This pattern implies that depletion of gas-phase coolants is likely to be the most 



important aspect of dust in terms of its influence on the observed spectrum (see also Matins 
1986); this statement is supported by comparisons of how forbidden line strengths change when 
the three effects of dust are added in isolation, rather than in succession. While depletion of 
optically emitting elements can produce a reduction of line emission at low Z, the result is more 
typically an increase in optical line emission as the cooling efficiency of other transitions decreases 
and equilibrium T e increases. Infrared fine-structure transitions are particularly important for 
establishing this behavior, since they scale in strength with abundance but are relatively insensitive 



to T e . Henry| (1993) has emphasized the importance of the [Sill] 34.8 fim line and infrared 



[Fell] transitions for cooling in an undepleted nebula; in real nebulae, these elements are likely 
to be subject to significant depletion onto grains, producing an increase in T e that can boost the 
strength of optical lines. 

Selective absorption of ionizing radiation by dust additionally modifies the ionization and 
thermal structure of a nebula, leading to changes in its optical spectrum. The decrease in grain 
opacity at photon energies above the Lyman limit results in a hardening of the ionizing radiation 
field as it is attenuated by grains. The resultant boost in average photon energy produces 
greater heating per ionization, with a corresponding increase in T e and optical line emission. The 
amplitude of this effect in GEHRs is likely to be small, based on the results in Figures 1 and 2 as 
well as analytic arguments. The optical depth of dust will scale with the total column density of 
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gas within the Stromgren radius, Njj = nnRs, and the dust-to-gas ratio, assumed here to scale 
with Z. From the definition of U, we can express 

Nh = U—— ~ 3 x 10 20 (JL) cm- 2 , (1) 

where the case B hydrogen recombination coefficient as for T e = 7500 K (approximately correct 
for Z = Zq) and n e = 100 cm -3 is taken from Hummer Storey) (1987). For a dust-to-gas ratio 



representative of the Galactic ISM, the maximum EUV absorption optical depth will occur at 
approximately 17 eV, with teuv ~ (2 x 10~ cm 2 )iVjj ( Martin Rouleau 1989). For a Stromgren 
sphere, the maximum EUV optical depth with thus be 

teuv ~ 0.6 (-^] , (2) 



V10-3; Vz 



where U is defined in the dust-free limit. The expected teuv m GEHRs is thus significant for 
typical parameters, and will be higher if the nebulae are relatively dusty and/or U is elevated 
above typical values. Substantial modification of the ionizing spectral energy distribution is more 
likely in Galactic ultracompact and blister HII regions (e.g., Aannestacf 1989), which may be 



described by U values much greater than those typical of GEHRs. When dust opacity is large, R s 
and teuv are obviously coupled further by the diminution of ionizing flux seen by atomic material, 



a point considered in more detail in §4.3. 



Figures 1 and 2 also demonstrate that grain heating boosts the predicted optical emission even 
further. Photoelectrons emitted by grains are expected to be more energetic and hence produce 
greater heating than photoelectrons emitted by the constituent atoms in isolation, because the 



ionization potential, or work function, is reduced for matter in a solid phase (e.g., see Draine 1989 
for a review). The treatment of grains employed by Cloudy assumes a minimum work function of 
8 eV for a neutral grain, with an increase in threshold energy as the grain becomes more strongly 
depleted of electrons and charged. The heating rate scales directly with the density of grains and, 
for a weak grain potential, the flux of radiation above the threshold energy. In a dust-free nebula, 
heating tends to be dominated by photoelectrons from recombined hydrogen and hence scales as 
ra 2 . For an HII region with fixed density and dust-to-gas ratio, the significance of grain heating 
relative to hydrogen bound-free heating consequently scales with the flux <p of energetic radiation 
(or more generally, with ionization parameter U, defined locally as a ratio of ionizing photon to 
electron densities) . The relative dependence of grain heating on introduces a radial variation in 
the resulting thermal effects for a spherical nebula, with greater heating occurring at small radii 
(see § |5,2.|) . Grain cooling by electron attachment is also likely to be significant, and the relative 
importance of heating and cooling by grains is a more complicated function of the ionizing spectral 



energy distribution, grain potential, and T e (Baldwin et al. 1991, Appendix C.14). 



These initial results suggest that the collective effects of dust may indeed account for the 
absence of observed pure Balmer-line nebulae. While [O III] emission may become unobservable 
in low excitation GEHRs, [O II] emission is invariably detected, as noted previously, with [O II] 
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A3727/H/3 ^ 0.3. This minimum is consistent in the present calculations with Z as high as 
2 — 3x Z . Observed strengths of [Nil] and [S II] are similarly consistent in general terms with 
metallicities as high as this bound. 



4. Verisimilitude 



4.1. Grain Assumptions 



The calculations presented here make a number of assumptions of uncertain validity for 
GEHRs. The greatest uncertainty in this discussion is the detailed nature of grains in these 
nebulae. From studies of the Galactic ISM, it is well established that depletion rates are a function 
of environment ( |Cowie Sz Songaiia 1986 and references therein). Given the substantial output of 
radiative and kinetic energy by young stars, grains located near these sources might be expected 



to be modified and/or partially destroyed. Strom et al.| (1974) presented a comparison of infrared 
emission from Galactic and extragalactic HII regions that has been interpreted as evidence for 
grain destruction in the GEHRs. Galactic sources exhibit a correlation between thermal radio 



continuum emission and infrared luminosity that Strom et al. employed to predict infrared 



fluxes from a group of more luminous GEHRs, using H/3 emission in place of radio continuum 
emission to trace Q(H). Observed 10 — 20 /mi fluxes for the extragalactic sources fall short of the 
predicted fluxes derived in this manner. [Strom et al. discussed several possible explanations for 
this discrepancy, including a deficiency of dust in the luminous extragalactic sources. 



The difference noted by [Strom et al. between measured IR fluxes and the extrapolated 
predictions may be a manifestation of a more general pattern of disparity between high- and 
low-luminosity sources. As discussed by Habing fc Israej (1979) and Kennicutt] (1984), HII regions 
show systematic trends with luminosity such that lower luminosity sources are denser and more 
centrally concentrated than their giant counterparts. Both tendencies might be expected to result 
in higher U for the Galactic sources, and consequently greater heating of dust grains and resultant 
IR emission, consistent with the observational trend found by Strom et al.. Alternatively, studies 



of the Galactic ISM demonstrate that depletion and hence the dust-to-gas ratio correlate with 
density ( Mathis| 1990), which would be consistent with a higher grain content in dense Galactic 
HII regions when compared with more diffuse GEHRs. 

Recent clues to the grain content of GEHRs are provided by analyses of silicon and iron 



abundances in HII regions and compact blue galaxies by Garnett et al. (1995) and Thuan (1995, 
private communication). In both studies the results suggest that the depletion rates of these 
elements are lower in GEHRs than in the average ISM. These findings support a picture in which 
a fraction of the grains in nebulae are destroyed, such that the dust-to-gas ratio in GEHRs is 
reduced somewhat from a simple extrapolation based on the Galactic ISM. If these results are 
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characteristic of GEHRs, the nebular predictions described here will exaggerate to some degree 
the actual effects of dust in these sources. 

The present calculations also assume a specific grain size distribution, which is again known 
to be a function of environment. The EUV grain opacity and rate of grain heating are very 
sensitive to the small grain component; if the small grain population is reduced compared to the 
ISM average, as inferred in the Orion nebula and other dense environs ( jMathis 1990 and references 



therein), the effects of dust opacity and heating will be correspondingly diminished. If 30 Doradus 
can be regarded as representative of GEHRs, however, the small grain component may actually be 
enhanced relative to grains in the Galactic ISM (Fitzpatrick 1986). The elemental content of dust 
in GEHRs is a further major uncertainty. The relative depletion rates are poorly constrained, and 
more generally, the total relative abundances may differ from a simple scaling of solar proportions. 
In light of these many fundamental ambiguities, the present calculations can only be considered 
illustrative of the effects of dust in HII regions, rather than providing definitive predictions. 



4.2. Line Ratios 



One means of assessing whether the nebular calculations described here bear any applicability 
to real GEHRs is by comparison of predicted and observed line ratios for these systems. Figure 3 
shows line ratio plots selected for ease of comparison with those given in Figure 8 of the extensive 
study of GEHRs by McCall et al. (1985); the plotted points represent observational measurements 



from McCall et al. . The plotted curves show the prediction results for a dust-free nebula (solid 
curves) and for a depleted nebula with full grain physics (dotted curves) , as a function of changing 
Z, for the two choices of T*. The behavior of the oxygen lines suggests that a higher T* is required 



at lower Z, in agreement with previous work (e.g., McCall et al.| 1985). 



Figure 3 shows that in the high metallicity regime (Z ^ Z Q ), the calculations for the 
undepleted nebula tend to underpredict emission in [Nil] and [S II] , while the flux in these 
lines is excessive in the calculations with grains. This pattern is more readily seen in Figure 
4, which plots [OlIl]/H/3 as a function of [Nil] /Ha and [Sll]/Ha for the same data sets. A 
simple explanation for this behavior might be that the assumed grain component is excessive, and 
the actual depletion and grain effects are somewhat reduced. Other explanations are possible, 
however. The strength of [S II] emission is particularly sensitive to the value of U (e.g., (Bkillman 



1989), and additional calculations indicate that the predicted [Sll]/Ha ratio for the dusty nebula 
can be brought into agreement with the observations by increasing U by a factor of ~ 2 — 3. This 
remedy is less successful for the [N II] line strength, however, and fine-tuning of other parameters 
such as N abundance is likely to be necessary to bring predictions for [Nil] into agreement with 
the observations. Decreasing U in the dust-free case does not significantly increase [S II] , and leads 
to substantial decreases in the strength of the other optical forbidden lines considered here. 
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The predictions in Figure 4 permit an interesting comparison with similar plots in Figures 7 
and 8 of Kennicutt, Keel, & Blaha (1989) for observations of HII regions in the nuclei of galaxies. 



Kennicutt et al.| (1989) noted that the nuclear sources show a systematic trend of greater [Nil] 



and [S II] line strength at a given value of [OlIl]/H/3 than is seen in extranuclear GEHRs. Their 
observational results for the nuclear sources are, in fact, in relatively good agreement with the 
predictions for dusty nebulae presented here, for Z up to ~ 2Z Q . |Kennicutt et al, suggested that 
the trend in the nuclear H II regions might stem from the presence of weak active nuclei that are 
spatially unresolved from normal circumnuclear star- forming regions. The present results suggest 
a possible alternative explanation whereby the enhanced [N II] and [S II] emission results from an 
elevated dust content in nuclear nebulae relative to disk HII regions. The nucleus corresponds 
to the center of a galaxy's potential well, and the interstellar pressure in such an environment is 
expected to exceed typical disk ISM pressures by up to 2 orders of magnitude or more (e.g., Spergel 
& Blitz 1992; Heifer & Blitz 1993). Resulting enhancements in density may foster enhanced grain 
formation and growth, consistent with local correlations between density and depletion, thus 
providing a physical motivation for why nuclear nebulae could be unusually dusty. 



4.3. Observable Extinction 



If grains are important in modifying the overall behavior of an H II region, we might expect 
this fact to be signaled observationally by signatures of reddening or extinction. The extinction 
optical depth produced by dust associated with the HII scales as in equation 2, but with optical 
depth at V (iy) reduced by a factor of ~ 5 from teuv ( [Martin fc Rouleau] 1989), such that Ty 
will be only a few tenths as the dust column becomes optically thick at ionizing energies. 

If the dust optical depth is increased beyond this level, the fact that the grain EUV and 
optical opacities are coupled implies an additional feedback mechanism that will limit optical 
indications of extinction. Significant teuv reduces Nh) but achieving a large teuv requires a 
substantial Nh for typical dust-to-gas ratios. Large EUV opacity and hence significant Ty will 
thus result only from extreme nebular conditions. For example, if internal dust is sufficient to 
produce Ty w 1 and hence teuv ~ 5, then Nh (which is proportional to [Q(H)ne 2 ] lj/3 in the 
dust-free case) must be ~ 3 x 10 20 for Z = Z Q (cf. eqn. 1). But Q(H) seen by the gas will 
be reduced as teuv increases and a growing fraction of the ionizing radiation is absorbed by 
grains; teuv ~ 5 will result in a reduction of the nebular radius by a factor of ~ 2 relative to the 
dust-free R s , with ~ 90% of the ionizing photons absorbed by dust for a standard dust-to-gas 
ratio ( Petrosian, Silk, Field 1972; see also Mathisj 1971 and Aannestad 1989). In order to obtain 



the Nh producing « 1, the product [Q(i?)ne 2 ] 1 / 3 and hence U corresponding to the dust-free 
case must be increased by approximately an order of magnitude. 



Conditions in which dust so strongly modifies the nebula are probably only present in systems 
that are clearly distinct from classical GEHRs. One such class of objects is that of ultracompact 
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HII regions (e.g., Wood fc Churchwell| 1989), which are far denser than typical GEHRs and show 



observational evidence for substantial dust absorption of ionizing radiation. Similar conditions 
may also obtain in the dense central regions of starburst galaxies. The high n e and large U that 
describe these environments are inconsistent with optical line ratios for GEHRs. 

Substantial extinction may still be present with less dramatic effects on the nebula if the 
ionization is produced by distributed sources rather than a central compact star cluster. An 
idealized picture illustrating this point would be a collection of Stromgren spheres aligned along 
the line of sight; ry through a single sphere may be small, but the collective effect of multiple 
foreground systems could produce significant extinction of emission produced in the more distant 
spheres. 

Observations suggest, however, that internal optical extinction in GEHRs is indeed fairly 
small. Studies of the stellar population in 30 Doradus ( Malamuth fc Heap| 1994) and in NGC 3603 



(Moffat, Drissen, & Shara 1994), for example, demonstrate that the massive stellar component that 



dominates production of ionizing radiation is relatively concentrated, suggesting that idealization 
of the nebula as a single Stromgren sphere is unlikely to be grossly in error. Empirical estimates of 
Ay for GEHRs are also rarely more than Ay ~ 1 ( Kennicutt) 1984). While extinction values based 



on reddening of the Balmer lines alone may underestimate the total extinction, due to surface 
effects ( |Mathis| 1970), radio measurements of GEHRs suggest that the optical estimates of total 
Ay are not in error by large factors ( |Israel &: Kennicutt 1980; C apian &: Deharveng 1986; van der 



Hulst et al. 1988; see also discussion in |Shields| 1990). 



Although the total optical extinction in GEHRs is thus only modest, the effect of dust on the 
nebula can still be significant if EUV opacity is much larger than opacity at optical wavelengths. 
It also must be remembered that attenuation of the ionizing radiation field is not directly traced 
by the comparison between Balmer line and radio fluxes. 



5. Further Implications 



5.1. Global Metallicity 



Existing abundance analyses of HII regions have relied heavily on the R23 parameter [= 
([Oil] A3727 + [OIII] AA4959,5007)/H/3] as an indicator of the abundance of oxygen, which in 
turn is often taken as a proxy for the total metallicity flPagel et al.| 1979). This approach can be 
justified by the dominant role of oxygen as a coolant in typical nebulae; the calibration of R23 
as a function of Z is nonetheless dependent at some level on the abundances of other coolants, 
however. An important question is thus whether this calibration is significantly modified by the 
presence of grains and associated depletion. Another way of expressing this question is to ask 
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whether the gas-phase oxygen abundance derived from R23 can be translated into a total oxygen 
abundance (or Z) through simply scaling by the oxygen depletion factor. 

To illustrate the effects of dust on the use of -R23 as an abundance indicator, Figure 5 displays 
R23 as a function of oxygen abundance in the gas phase, for the calculations described in § j2l By 
comparing the curves for a given T± that correspond to no dust and full dust effects, we can see 
that relatively small changes in R23 result from the presence of grains and associated depletion 
at low Z. The shift in predicted R23 due to grains at a given gas-phase abundance increases 
significantly at higher Z; at the same time, however, R23 is becoming a steeper function of Z, 
such that a large shift in predicted R23 propagates into a small difference in Z. In terms of the 
observationally relevant comparison, for a measured value of R23, the uncertainty in Z due to dust 
effects is likely to be ^ 0.2 dex for a gas-phase oxygen abundance of up to ~ 3 times the solar 
value. We conclude that the effects of grains are unlikely to introduce large errors in measurements 



of gas-phase oxygen abundance based on R23 . This conclusion differs somewhat from that of Henry 
(1993). The calibration of R23 at high Z is rather sensitive to n e , however, due to suppression of 
fine-structure cooling ( Oey fc Kennicutt| 1993), as well as U, and the cumulative effects of these 



uncertainties cast doubt on whether a reliable calibration of R23 versus Z is possible at Z > Z 
(see also |Mathis| 1986). 



5.2. Temperature Gradients 



Temperature gradients in dust-free H II regions are determined by the efficiency of cooling 
mechanisms and the average heating per ionization as a function of radius. At very low Z, 
collisionally excited hydrogen Lya plays an important role in cooling. The neutral fraction of 
hydrogen that can give rise to such emission increases monotonically with radius, so that the 
cooling efficiency increases with r, and T e tends to be a decreasing function of radius. At higher Z, 
collisionally excited lines of heavy elements become important for cooling, with infrared transitions 
playing an increasing role; important ions with ground-state fine structure giving rise to infrared 
emission include doubly ionized O, S, and N, which boosts the cooling efficiency in the inner 
nebula where these species are abundant^ (|Stasihska| 1980). In addition, the frequency dependence 
of hydrogen and helium bound-free opacity guarantees that the average energy and resultant 
heating per ionizing photon increases with increasing radius. The combined effects are expected to 
produce a gradient dT e /dr > in high metallicity HII regions^. Infrared transitions also play an 
important role in cooling the lower-ionization, outer portions of a nebula, particularly by emission 
of [Sill] 34.8 jim., lines of [Fell] (e.g., see Nussbaumer k, Storey| 1988), and [Cll] 158 /xm. 



2 The resulting lines are [O III] 51.8 and 88.4 ^m, [S III] 18.7 and 33.5 fim, and [NIII] 57.3 /xm 

3 At very low Z the gradient in heating stemming from bound-free opacity is overwhelmed by the role of Lya 
cooling in determining the equilibrium T e , which consequently decreases with r 
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Grains and associated depletion may have an important role in modifying the thermal 
structure of GEHRs beyond a simple global increase in T e , when compared with the dust-free 
case. Examination of the elemental depletion factors given in Table 1 indicates that depletion 
has a much stronger influence on elements dominating fine-structure cooling in the low-ionization, 
outer portions of a nebula (Si, Fe, C) than on those important in the inner nebula (O, S, N). 
Depletion will consequently tend to steepen temperature gradients from the undepleted case, by 
preferentially elevating the temperature of the outer nebula. This steepening will be mitigated, 
however, through heating of the inner nebula by grains formed from the depleted elements. As 
noted in §^], the rate of grain heating will tend to scale locally with the flux of ionizing radiation. 
Radial dilution and attenuation of the radiation produced by a central star or stars will thus lead 
to a radial falloff in grain heating within a nebula, assuming a uniform distribution of dust. 

These consequences can be seen quantitatively in the middle panels of Figures 6 and 7, which 
show T e profiles for the calculations described in §^], for the case of Z = Z . The results of adding 
the effects of dust incrementally are illustrated by the successive curves. For comparison, the lower 
panels show the characteristic ionization state of the gas as a function of radius, as described by 
the ionization fraction of + and ++ . Grain heating and cooling as a function of r is illustrated 
in the top panels. Grain heating may exceed 70% of the total heating rate in the inner nebula, 
and decreases monotonically with increasing radius through most of the HII region. Near R s , the 
fractional heating may again become large; in this partially ionized zone, hydrogen bound-free 
heating is rapidly diminishing, but photons with energy 8 ^ hv < 13.6 eV continue to produce 
significant grain photoelectrons. Fractional cooling by dust is lower than fractional heating at 
all radii, with the possible exception of the outermost portions of the Strdmgren sphere. These 
results show a general similarity to calculations with a blister geometry described by [Baldwin~et 
|I| (1991). 

The radial dependences of depletion and grain heating propagate into the predicted fluxes of 
optical lines, as can be seen by a review of Figures 1 and 2. For both choices of T*, emission in 
[Oil], [S II], and [Nil] (which represent relatively weakly ionized species) show a strong response 
to depletion, while the response in [O III] is less dramatic. In contrast, [O III] shows a much greater 
sensitivity to grain heating than do [Oil], [S II] , or [Nil], for the nebula powered by the cooler 
star. The differences in the behavior of [O III] between the two choices of T* reflect the relative 
extent of the ++ zone (cf. Figs. 6 and 7). The cooler star produces a more concentrated region 
of ++ , which is consequently more strongly affected by grain heating; the ++ zone produced by 
the hotter star overlaps more broadly with zones of low-ionization species, introducing a greater 
sensitivity to depletion. 

Grain heating might be expected to produce very strong enhancements of high-ionization lines 
that form close to the ionizing stars. In the present calculations, however, an easily observable 
signature of this effect is not predicted, because such lines are relatively weak, even when dust 
heating is taken into account. A tradeoff results in that a strong dust signature may be more 
apparent when T* is relatively low, due to the fact that the emitting volume of a high-ionization 
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species is then highly concentrated to the region of strong heating by grains; but when this volume 
is small, emission is likely to be weak except for a tracer with high abundance, such as oxygen. 
An illustration is provided in the case of = 38,000 K by [Nelll] A3869/H/3, which increases 
by a factor of ~ 2 at Z = Zq and by an order of magnitude at Z = 2Z Q when grain heating 
and cooling is introduced to the calculation; but invariably the [Nelll] line flux is less than 1% 
of that of H/3. If T+ is as high as 45,000 K, the [Nelll] line becomes more readily observable 
([NeIIl]/H/3 = 0.1 — 0.2), but shows only a weak sensitivity to dust thermal effects since the Ne ++ 
volume extends over a larger fraction of the nebula, including regions where grain photoelectric 
heating is relatively unimportant. 

Signatures of radially dependent grain heating may be more readily detected in planetary 
nebulae, which can span a broader range of ionization than is typically found in HII regions. 
Keyes, Aller, fc Feibelman| (1990) performed a detailed analysis of the spectrum of the planetary 



nebula NGC 7027, and noted a persistent difficulty in accounting for the observed (strong) 
emission in [NeV] AA3346, 3426 after adjusting model nebular parameters so as to reproduce 
correctly the strength of He II lines. Since the He II features are recombination lines while [Ne V] 
is collisionally excited, the unexpected strength of the latter emission might be explained by a 
thermal mechanism such as grain heating in the inner nebula, if a fraction of the dust survives in 
this environment. Test calculations with model parameters similar to those employed by Keyes 
et al. are consistent with this interpretation. Detailed study of other sources is needed in order 
to ascertain whether grain heating is often an important influence on the emergent spectrum of 
planetary nebulae. 

Abundance analyses based on HII region spectra are potentially sensitive to the details of 
nebular temperature gradients. The derivation of an ionic abundance from a collisionally excited 
emission line requires the use of a T e value appropriate to that ion. Nonuniform temperature 
structure necessitates the use of different characteristic temperatures for different ions; in addition, 
temperature-sensitive (auroral/nebular) line ratios of a single ion can differ from ion- weighted 
temperatures that are physically appropriate for extracting abundance information from line 
strengths. Radial gradients are a specific case of temperature fluctuations that contribute lingering 



uncertainty in the calibration of abundances based on nebular studies (e.g., Peimbert 1967; 



Peimbert, Storey, Sz Torres-Peimbert] 1993). Detailed aspects of these problems have been reviewed 



by |Garnett (1992), who derived theoretical relations between ion- weighted temperatures for 



various species, as well as between line-ratio and ion-weighted temperatures. Garnett's models do 
not include the consequences of grains, although he notes their possible significance for perturbing 
nebular temperature structure. 

The behavior of ion-weighted temperatures in the present calculations suggest some possible 
modifications to previous theoretical work. Garnett's models included depletion of Si (which he 
excluded completely) and Fe, and our calculations with only depletion are consistent with his 
results. In the preceding discussion, it was noted that depletion selectively boosts T e in the outer 
nebula, while grain heating is most important in the inner nebula. As a result, models that include 
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depletion but neglect the thermal effects of the corresponding grain population may overestimate 
the steepness of the temperature gradient in an enriched nebula. These expectations are borne out 
in part in the present calculations, as can be seen by examination of Figures 8 and 9, which plot 
the ion-weighted temperatures T(0 + ) and T(S ++ ) as functions of T(0 ++ ). These temperatures 
are calculated according to equation (1) in |Garnett| (1992); the plotted points span Z = 0.1 — 3 
Zq. Also shown in these figures are the trends found by Garnett (and confirmed by us) for a 
nebula with depletion alone, i.e., his equations (2) and (3), as well as diagonal lines of uniform T e . 

At high Z, the calculations with depletion and grain effects show a more nearly isothermal 
behavior than predictions restricted to depletion, for the case of the cooler T*. For comparison, the 
predictions for a dust-free and undepleted nebula are also shown; interestingly, these predictions 
are in better agreement with the calculations with full dust effects than with those taking only 
depletion into account. This behavior reflects a rough tendency for depletion and grain heating 
to balance each other in terms of their net effect on dT e /dr. For the hotter T*, T(S ++ ) at a 
given T(0 ++ ) for the dusty nebula is still somewhat overpredicted by the relation with depletion 
only, while T(0 + ) versus T(0 ++ ) is consistent with the depletion-only results. For T* = 45,000 
K, S ++ occupies an intermediate zone in radius that experiences a smaller relative increase in 
T e than either of the inner or outer regions weighted respectively by ++ or + ; the overall 
gradient in T e remains positive for the most part, however, such that the relation between T(0 ++ ) 
and T(0 + ) does not change appreciably from the dust-free or depletion-only case. Calculations 
with n e increased to 100 and 1000 cm -3 produce a general increase in nebular temperatures, 
due to partial collisional suppression of fine-structure cooling, but the ion-weighted temperature 
behavior is qualitatively the same as seen in Figures 8 and 9. These results also generally show 
a weak sensitivity to U. Increasing e such that log U = —2.5 causes the dust-free calculations 
for T* = 38, 000 K to approach Garnett's relations, since the inner zone of high-ionization that 
cools efficiently expands, resulting in a steeper T e gradient; the results for the dusty nebula 
are qualitatively unchanged, however. We find reasonable agreement in the present calculations 
with relations listed by |Garnett| for other ion-weighted temperatures [i.e., T(Ar ++ ) = T(S ++ ); 
T(N+) = T(0+); T(N++) = T(0 ++ ); T(Ne++) « T(0++); T(S+) « T(0+)]. 

If enriched GEHRs contain significant quantities of dust, the present results suggest that these 
nebulae may be treated as isothermal systems in terms of the relative behavior of ion-weighted 
temperatures, to good approximation. While this statement is strictly true only for the calculations 
presented here with the cooler T*, empirical evidence suggests that the characteristic T* of a 
nebula decreases with increasing metallicity, such that a relatively soft radiation field is, in fact, 
appropriate for enriched systems with T e £ 7000 K. If dust is efficiently destroyed in HII regions, 
an isothermal treatment may still be appropriate, since the effects of depletion as well as grain 
heating would be vitiated^. 



4 A relevant technical point is that Cloudy employs a modified on-the-spot (OTS) approximation for transfer of 
the hydrogen and helium diffuse radiation field, fn some cases OTS can lead to an estimated temperature gradient 
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We note that nebular models with depletion alone may nonetheless produce accurate 
predictions if grains are preferentially destroyed at small radii. In this situation the effects of grain 
heating that predominate at small r would be removed, while depletion effects that selectively 
influence the outer nebula would still be relevant. Evidence for reduced depletion in the inner, 



high-ionization regions of planetary nebulae has been discussed recently by Kingdon, Ferland, & 
Feibelmanl (1995). 



Modifications to relations between line-ratio-inferred T e and ion-weighted temperatures will 
also result from dust effects. In the present calculations the largest discrepancies between the 
two temperature values occur at ion-weighted values of T e £ 6000 K; in this domain, auroral 
transitions employed for direct estimation of T e are generally too weak to be observable. Additional 
perturbations in temperature estimates from line ratios (and also ion-weighted T e values) may 
result from more complicated spatial distributions of dust, as noted by Peimbert et al. (1993). 

As a final cautionary note, we remind the reader that the present calculations (and many 
others in the literature) assume a constant density, independent of r. While a detailed inquiry 
into the effects of nonuniform density is beyond the scope of the present work, line fluxes and 
temperature structure can certainly deviate in important ways from the predictions given here if 
large gradients in density exist within a nebula. Observations suggest that n e is often a decreasing 
function of r in GEHRs (e.g., Kennicutt| 1984). For HII regions with gradients of this form, the 



temperature of the inner nebula may be boosted and positive gradients in T e again diminished if 
n e at small radii is sufficient to partially suppress cooling via infrared fine-structure transitions 
with low critical densities. In general terms, however, the physical effects of the generic grains 
discussed in this paper are independent of the density structure of the nebula in which the dust 
resides. 



6. Conclusions 



Embedded dust can modify the optical spectrum of GEHRs, through a combination of effects 
that mostly act to increase T e . The largest perturbations are generally the result of depletion 
of coolants from the gas phase, although grain opacity and heating may be significant in some 
cases. The net result is usually an increase in optical forbidden-line strengths, particularly at high 
metallicity. Dust may lead to nonnegligible modification of nebular behavior without introducing 
strong signatures of reddening or extinction. 



that is too steep compared to an outward-only approximation (which was employed by Garnett 
treat ment of radiativ e transfer (i.e., dT e /dr is too large; see, for example, discussion in [Barazin 



Ferland et al 



1992) or a rigorous 



1977; McCall et al 



1985; 

(positive) dT e /dr that is too shallow. 



1995), which suggests that our results are conservative in that we are unlikely to predict a 
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The influence of grains provides a simple way of explaining the observational lack of H II 
regions emitting a pure Balmer-line spectrum at optical wavelengths, which is expected on 
theoretical grounds for dust-free nebulae with metallicity only slightly greater than Z Q . At high 
Z, predicted enhancements of [S II] and [N II] in particular show a resemblance to the empirical 
behavior of H II regions in galaxy nuclei. If these environments are particularly dusty, due perhaps 
to characteristically elevated densities, the influence of dust provides an alternative mechanism to 
weak active nuclei for explaining the anomalous line ratios of these sources. 

Depletion and heating by photoelectrons emitted by grains in H II regions may introduce 
significant modifications to nebular temperature structure. Depletion selectively influences 
elements that are important for cooling the outer portions of nebulae, while grain heating is 
most effective at small radii. Modifications to temperature gradients within GEHRs by dust may 
necessitate changes in assumptions employed in nebular abundance analyses, specifically in the 
relationships between ion-weighted temperatures for different species. For chemically enriched 
GEHRs, the present results suggest that these nebulae can be treated as nearly isothermal systems. 

Inclusion of grain physics is unlikely to result in large changes of gas-phase heavy-element 
abundances derived for low-Z objects and GEHRs ionized by relatively hot stars. In general, 
dust-related phenomena will contribute an uncertainty of ~ 0.2 dex or less in gas-phase oxygen 
abundances derived from the -R23 parameter. Modifications to nebular temperature structure by 
dust and associated depletion may be significant in HII regions described by high Z and low T*, 
and should be considered when deriving relative abundances for heavy elements. 
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Telescope Science Institute, which is operated by the Association of Universities for Research 
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Informative discussions with John Black, Don Garnett, Sally Oey, and Evan Skillman are 
gratefully acknowledged. We also thank the anonymous referee for helpful comments on the 
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TABLE 1 



Relative Abundances 21 at Z = Z Q 



Element 


Total 


Depleted 


He 


10.99 


10.99 


C 


8.56 


8.22 


N 


8.05 


7.95 





8.93 


8.73 


Ne 


8.09 


8.09 


Na 


6.30 


5.50 


Mg 


7.58 


7.15 


Al 


6.47 


5.53 


Si 


7.55 


6.55 


S 


7.21 


7.21 


Ar 


6.56 


6.56 


Ca 


6.36 


3.66 


Fe 


7.67 


5.97 


Ni 


6.25 


4.26 



a Abundances are logarithmic + 12, by number relative to hydrogen. 
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Figure Captions 

Fig. 1 — Line ratios plotted as a function of total metallicity, with dust effects added 
incrementally, for T* = 38, 000 K. Solid curves: no depletion or grains; dotted: depletion only; 
short dashed: depletion + grain opacity; long dashed: depletion, with grain opacity, heating, and 
cooling. The lines [Oil] = AA3726,3729; [OlII] = AA4959, 5007; [S II] = AA6716,6731; [Nil] 
= AA6548, 6583. Approximate lower bounds to observed ratios in HII regions are indicated by the 
dotted horizontal lines. 

Fig. 2 — Same as Figure 1, for T* = 45, 000 K. 

Fig. 3 — Loci of line ratios resulting from variation in Z (= 0.1 — 5 Zq), for dust-free (solid 
curves) and dusty (dotted curves) nebulae. Values of Z are indicated by the corresponding filled 



points. Crosses represent measured, dereddened values for GEHRs taken from McCall et al. 
(1985). 

Fig. 4 — Same as Figure 3, with axes chosen for comparison with results for nuclear and disk 
GEHRs studied by [Kennicutt et al.| (1989). A general discussion of the loci of emission-line nuclei 
and other objects in line-ratio diagrams can be found in Veilleux &; Osterbrock| (1987). 



Fig. 5 — Values of R23 as a function of gas-phase oxygen abundance. Results are shown for 
nebular cases that are dust-free (solid curves) and dusty (dotted curves), with oxygen abundance 
normalized to log (O/H)© = —3.07 by number. The upper two curves assume T* = 45,000 K, 
while the lower pair assumes T± = 38, 000 K. The abundances of other heavy elements scale with 
oxygen as described in §2. 

Fig. 6 — Radial behavior of thermal and ionization properties for T* = 38, 000 K and Z = 
Z . Central panel: T e as a function of radius, with depletion and grain effects added incrementally. 
Top panel: heating and cooling contributed by grains as a fraction of the total. Bottom panel: 
Ionization fraction of + and ++ ; solid curves represent a dust-free nebula, while the dashed 
curves show the result after adding full dust physics. 

Fig. 7 — Same as Figure 6, for T* = 45, 000 K. 

Fig. 8 — T(0 + ) versus T(0 ++ ) for the dust-free and dusty calculations, where the 
temperatures represent ion- weighted averages. The dashed line represents the theoretical relation 
obtained by Garnett] (1992) for depleted nebulae, while the dotted line represents the isothermal 
case. 

Fig. 9 — T(S ++ ) versus T(0 ++ ) for the dust-free and dusty calculations. The dashed line 
represents the theoretical relation obtained by |Garnett| (1992) for depleted nebulae, while the 
dotted line represents the isothermal case. 




»H/[ll S] &0| 



»H/[ll N] 6o| 




»H/[ll S] &0| 



»H/[ll N] 6o| 




R (cm) 




2x10 4x10 
R (cm) 



6x10 



1.5 — 



o 



1 — 



0.5 — 



o 



1 — 



0.5 



T. = 38,000 K 
a no dust 
a dust 



.A' 



,4 



' A 



1.5 — T, = 45,000 K 



s± / 

A 



& / 



0.5 1 1.5 

T(0 +2 ) / (10 4 K) 



1.5 — 



o 



1 — 



CM 
+ 

CO 



0.5 - 



o 



1 — 



CM 
+ 

CO 



0.5 — 



T. = 38,000 K 
a no dust 
a dust 



s 



s 



S' 



' A/ 
A- 



1.5 — T, = 45,000 K 



s 

s 



.'A'' 



A.- 



o 



0.5 1 1.5 

T(0 +2 ) / (10 4 K) 



2 



